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ABSTRACT 



We present measurements of the Fundamental Plane (FP) parameters (the effec- 
tive radius, the mean effective surface brightness, and the central velocity dispersion) 
^ ' of six field elliptical galaxies at intermediate redshift. The imaging is taken from the 

I Medium Deep Survey of the Hubble Space Telescope, while the kinematical data are 

CN ■ obtained from long-slit spectroscopy using the 3.6-m ESO telescope. The Fundamental 

' Plane appears well defined in the field even at redshift « 0.3. The data show a shift 

in the FP zero point with respect to the local relation, possibly indicating modest 
evolution, consistent with the result found for intermediate redshift cluster samples. 
The FP slopes derived for our field data, plus other cluster ellipticals at intermediate 
redshift taken from the literature, differ from the local ones, but are still consistent 
• with the interpretation of the FP as a result of homology, of the virial theorem and 

of the existence of a relation between luminosity and mass, L oc M^. We also derive 
Q . the surface brightness vs. effective radius relation for nine galaxies with redshift up to 

$-H ' -2 ~ 0.6, and data from the literature; the evolution that can be inferred is consistent 

, with what is found using the FP. 
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1 INTRODUCTION 

Elliptical galaxies in the local universe are known to pop- 
ulate only a two-dimensional manifold, known as the Fun- 
damental Plane (hereafter FP; Djorgovski & Davis 1987; 
Dressier et al. 1987), in the three-dimensional space defined 
by effective radius _Re, mean surface brightness within the 
effective radius SBe (SBe ~ — 2.5 log(7)e-|-const; I is the 
surface brightness in linear flux), and central velocity dis- 
persion a. The FP is described by: 

log i?e = a log (7 + /3 SBe -F- 7, (1) 

Re is in kpc (while re is in arcsecs), a in km s~^, SBe 
in mag arcsec"^ and Ho is hereafter assumed to be 50 /150 
km s~^Mpc~^. Typical values of the coefficients are for ex- 
ample a = 1.25, f} = 0.32 and 7 — —8.895 in Johnson B 
(Bender et al. 1998; hereafter B98). It is still unclear whether 
the thickness of the FP is intrinsic, i.e., due to scatter in the 



properties of elliptical galaxies, or is due to observational 
errors. 

The existence of the empirical scaling described by the 
FP has strong implications in terms of galactic evolution 
and formation theories. For example, the FP suggests that 
the distributions of dark and luminous matter are related, 
since a depends on the total gravitational potential, whereas 
SBe and -Re trace only the stellar mass. The FP can be ex- 
plained in terms of homology, of the virial theorem, and of 
the existence of a well-defined relation between luminosity 
and mass, L oc Af (Faber et al., 1987; see also van Al- 
bada, Bertin & Stiavelli 1995). This interpretation requires 
a relation between the derived values of a and /3, namely 
Q-10/3-f2 = (77 = 0.2 a//3). 

Measuring the FP parameters for ellipticals at interme- 
diate redshift z ~ 0.3-0.8 has only recently become feasi- 
ble (Pahre, Djorgovski & de Carvalho 1995; van Dokkum 
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& Franx 1996; Bender et al. 1996; Kelson et al. 1997; van et al. 1994). The target pool has been selected originally 



Dokkum et al. 1998). This opens up two important ques- 
tions which bear on the age, formation history, and inter- 
nal properties of elliptical galaxies: i) how far in the past 
does an FP-like relation apply? and ii) do its parameters 
evolve significantly with time? The existence of a tight FP- 
like relation at substantial look-back times would suggest 
the existence of a universal relation between mass and for- 
mation time for ellipticals, since, e.g., population differences 
would be amplified as the redshift increases. Existing mea- 
surements of the properties of ellipticals at intermediate red- 
shifts have concentrated so far on relatively rich c lusters of 
galaxies, where bright ellipticals are ea s ier to find ( Palireet 
al. 199|: pan Dokkum fc Franx 1991 iBender et al. igQ^ 



Ziegler _fc Bender 1997 ; Kelson et al. 1997 ; |yan Dokkum et 
al. 199^ and observations can be carried out efficiently with 
multi-object spectrographs. When comparing, e.g., the FP 
of intermediate redshift and zero redshift clusters one must 
take into account the evolution with redshift of the popu- 
lation of galaxies in clusters, and environmental effects. Al- 
though the ev olution of clusters of g alaxies is still somewhat 
controversial (Postman et al. 1996), the galaxy population 
in rich clusters is likely to become less uniform (in terms, 
e.g., of age and metallicity) with time, due to accretion of 
isolated galaxies and small groups. Moreover, the question of 
whether there are systematic differences between cluster and 
field ellipticals (the latter including galaxies in loos e groups 
or poo r clusters) is still o pen even at low redshift (le Car- 



valho & Djorgovski 1992). A study of the galaxy properties 



as a function of look-back time provides a sensitive probe 
of the possible evolutionary differences between cluster and 
field ellipticals. For these reasons, we have begun a study of 
the empirical scaling laws of field ellipticals at intermediate 
redshift. 

We have started out by defining a sample of ellipti- 
cals which is not biased in favour of rich clusters. Details 
on the selection process are given in Section ^ Sections ^ 
and ^ describe the data reduction and the error analysis. 
Measurements at non-zero redshift have to be carried out 
with photometric filters and effective apertures that differ 
from those used to derive the local scaling laws; Section ^ 
describes the conversion of our measurements to the stan- 
dard quantities. The main results are reported in Section ^ 
where our measurements are compared to local samples and 
intermediate redshift cluster samples. A summary is given 
in Section ^ 



2 SAMPLE SELECTION 

With this study we aim at addressing two questions: i) does 
the FP of field ellipticals evolve with redshift? and ii) do, 
at any given redshift, the field ellipticals have FP parame- 
ters different from those found in rich cluster environments? 
For this we need to select a primary sample of galaxies bi- 
ased against rich cluster membership at intermediate red- 
shift, where the project is feasible and the look-back time 
(« 4-5 Gyrs) is high enough that evolution can be notice- 
able. 

The targets used in this study have been chosen among 
a sample of random ellipticals found in the WFP C2 paral 
lei images collected by the Medium Deep Survey ( Griffith! 



from available HST images in the appropriate RA, 5 range 
according to the following criteria: 

(i) Apparent magnitude I < 20.5 

(ii) Morphology clearly defined as elliptical, with an ap- 
parent effective radius re > 0'.'5. 

(iii) V — I colour in agreement with the empiric relation 
for ellipticals in the redshift range 0.2 < z < 0.5 (1.1 < 
V-I < 1.7) 

Magnitudes and colours above are as defined by the 
Medium Deep Survey group (Griffiths et al. 1994); V and 
/ are the magnitudes in the WFPC2 filters F606W and 
F814W, computed in the WFPC2 Flight system (Holtzman 
et al. 1995), and the effective radius re is the value found 
from the two-dimensional image fitting carried out by the 
MDS. All photometric parameters have been rederived here 
with two different techniques (see Section ^ to ensure a 
self-consistent treatment. 

Individual targets were then selected from this pool of 
approximately 25 candidates on the basis of convenience and 
availability during each observing night; in the spirit of this 
exploratory study, we did not attempt further to achieve 
uniformity or completeness within our pool of candidates. 

Since we did not explicitly exclude members of poor 
clusters (no rich clusters were observed in these random 
fields), the ellipticals in our sample should be representative 
of a random, magnitude-selected sample of non-rich-cluster 
objects. In fact, four of the six objects for which we have 
obtained FP parameters happen to be members of a group 
(or a poor cluster), which was randomly observed in the 
available HST images. Future observations will enable us to 
achieve a better level of completeness and thus to better 
characterise our sample. 

Our colour selection criterion will bias against actively 
star-forming ellipticals, which ar e known to occ ur, albeit 
infrequently, in complete samples (Lilly et al. 1995). We plan 



to relax the blue cutoff in future observations. 



2.1 Target list 

The characteristics of the objects we observed are sum- 
marised in Table |^. The final sample used in the discussion in 
Section ^ is composed of the galaxies successfully fitted with 
the two different photometric techniques, isophotal profile fit 
and two-dimensional fit, and with measured kinematics. 

The photometry of C, L, M, N was of low quality due 
to companion galaxies, stars in the field of view and the 
faintness of the objects. The 2D fit was performed only on 
the galaxies with a high quality isophotal profile. We were 
able to determine high quality photometric parameters for G 
by subtracting a model of its companion H from the image 
before the 2D fit. 



3 PHOTOMETRY 

The images are taken fro m the Medium Deep Survey of the 
Hubble Space Telescope ( Griffiths et al. 1994 ) . For each field 
there are images from WFPC2 through filters F606W and 
F814W. Table |2| lists total exposure times and number of 
exposures. 
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Table 1. Galaxy identification. (1) Adopted names, (2) HST 
Medium Deep Survey field name, (3) WFPC2 chip, (4) ID within 
the field, (5) and (6) pixel positions of the center of the galaxy, 
(7) Sizes of the F606W and F814W images used in the 2D fits 
and displayed in Figures 2 and 3. The image used for I in F814W 
was 5"4 X 4'.'8. Dashes in Column 7 indicate that the 2D fit was 
not performed. 



galaxy 


field 


chip n 


X 


y 


size 


(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


A 


ut800 


4 


1 


324 


444 


18 X 18 


B 


ut800 


4 


2 


265 


741 


8" X 8" 


C 


ur610 


3 


1 


678 


691 




D 


ur610 


3 


2 


723 


642 


9'.' 5 X 6'.'4 


E 


u5405 


2 


1 


335 


324 


8" X 7" 


F 


u5405 


2 


2 


332 


221 


5'.'4 X 6'.'4 


G 


u5405 


2 


3 


192 


288 


10" X 10" 


H 


u5405 


2 


4 


194 


255 




I 


u5405 


2 


5 


72 


63 


5'.' 2 X 4'.'3 


L 


u5405 


3 


6 


254 


711 




M 


u5405 


3 


7 


266 


715 




N 


u5405 


3 


8 


747 


732 




O 


ustOO 


3 


1 


226 


523 


5" X 5" 


P 


ustOO 


3 


2 


336 


404 


5" X 5" 


Q 


ustOO 


3 


3 


443 


298 


5'.' 2 X 5'.'4 


R 


urzOO 


2 


1 


719 


242 




Table 2. Photometric data: number of exposures (nexp) and total 


exposure times (toxp). 










field 


filter 


nexp texp (s) 






u5405 


F606W 




1 


800 








F814W 




1 


800 






ur610 


F606W 




2 


1500 








F814W 




2 


1600 






urzOO 


F606W 




3 


5400 








F814W 




4 


8400 






ustOO 


F606W 




10 


16500 








F814W 




11 


23100 






ut800 


F606W 




2 


1430 








F814W 




3 


4430 







3.1 Reduction 

The basic reduction was done at the Space Telescope Science 
Institute using the standard pipeline with the best reference 
files available. As bias and dark files we used t he superbias 
and the superdark from the Hubble Deep Field (Williams et 
al. 1996) for their high signal-to-noise ratio. The cosmic ray 
removal for the fields with multiple images was done using 
the iraf/stsdas task CRREJ, with a cutoff at 5 standard 
deviations from the mean. Since only one exposure per filter 
was available for u5405, the original image of this field was 
compared with a smoothed version (using the MIDAS com- 
mand filter/cosmic). In all cases we checked by eye that 
the pixels removed in the area of the targets were only those 
affected by cosmic rays. 

For u5405, chip 2 of WFPC2 is aff'ected by a bad col- 
umn, very close to the center of galaxies E and F (see Fig- 
ure y) . We interpolated through it in the outer parts of the 
galaxies, but we chose not to modify the 4-5 pixel diameter 
central zone because it was too steep for any interpolation to 
be meaningful (see also Figure ^ . In order to test how much 



this bad column affected the determination of the photomet- 
ric parameters, we fitted the isophotal profile of E excluding 
the innermost part of the profile. We tried excluding a 2 pixel 
and a 4 pixel diameter circle; the effective radius and the ef- 
fective surface brightness change significantly (up to 20% 
in re), towards the value given by the 2D fit (smaller re 
and brighter surface brightness). These new values are still 
consistent with the average value within the internal errors 
(that are the largest of the sample, see Subsection i.4). Fur- 



thermore it should be noticed that the derived values for 
re and SBe are highly correlated and therefore the varia- 
tion of the combination that enters the FP, log re — /3SBe, 
is much smaller than the variation on the single coefficients 
(see Subsection 3^ and Table ^) . 

A constant sky contribution was subtracted from each 
image. The sky contribution used was the mean of the aver- 
ages of areas with no evidence for hght sources. This proce- 
dure does not take into account smooth and diffuse back- 
ground sources, which is one of the main sources of un- 
certainty when using curve of growth technique (see Sec- 
tion 



3.2 Isophotes and two-dimensional fits 

The photometric parameters needed for the FP are the ef- 
fective radius re, and the effective surface brightness SBe, 
measured as the best r^^* light profiles parameters. 

In order to obtain the best results and robust error 
estimates, the photometric parameters were derived using 
two independent techniques: fits to the isophotal luminosity 
profiles and two-dimensional fits to the images. Both proce- 
dures require a PSF to convolve the models wi th. We used 
synthetic PSFs calculated with Tiny Tim 4.0 ( [Krist 1994| ) 
with a 15 mas jitter. The quality of the synthetic PSF was 
checked by comparing the isophotal profiles of two stars and 
the profiles of two PSFs created in the same spot. No sig- 
nificant difference was noticed. As a double-check we fitted 
an isophotal profile using both a synthetic PSF profile and 
a real star profile: the differences were less than 1% in re 
and less than 0.02 mag in SBe. The photometric parame- 
ters of each galaxy were derived independently in the two 
bandpasses F814W and F606W. 



3.2.1 Isophotes 

Isophotes were fitted to the images with the center of the 
ellipse, the semiaxes a and b, and the position angle as 
free parameters; variations of ellipticity (e) and position an- 
gle (pa) with the "circularised" radius \/a& were allowed. 
Isophotes were derived using a version of the MIDAS com- 
mand fit/ell3 modified to better deal with steep gradients 
in the luminosity profile in the innermost pixels (M0ller, Sti- 
avelli & Zeilinger, 1995). The profiles obtained were fitted 
with an exponential law, an r^^'^ law, and a linear combina- 
tion of the two. In all cases, except for B, we found that the 
objects were best fitted by the r^'^^ law. A slightly better 
fit of the light profile of galaxy B was obtained by adding 
a small exponential component to the r^^* law, suggesting 
that the bump on the residual might be due to a fainter disk 
component (see Figures ^ and |^ . The fit was done usi ng th e 
least squares fitting software used by CaroUo et al. ( 1997 ). 
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Figure 1. Isophotal profiles. Crosses are the data points, solid lines represent the best-fit r^/* profile, convolved with the PSF. The 
upper data refer to F814W, the lower to F606W. The error bars (which are smaller than the crosses for the innermost part of the profile) 
take into account the sky-subtraction error. The bump on the B profile is possibly the signature of a disk component (see text, Section 



3.2.1 



and the two dimensional residual in Figure . The central surface brightness decline in E is likely to be caused by the bad column 



right at the center of the galaxy. 



For each galaxy we used a specific PSF calculated at its po- 
sition. In Figure |l| we show the best-fit r^'''* laws convolved 
with the PSFs, superimposed on the data points. 



3.2.2 Two-dimensional fits 

This technique consists in fitting a two-dimensional model of 
an elliptical galaxy, as seen on the CCD chip, to the WFPC2 



image. The model galaxy has an r^'''* luminosity profile, with 
fixed position angle and ellipticity (we used the value derived 
from the isophotal fit near the efi'ective radius) . Each fit has 
four free parameters: the efi'ective surface brightness, the 
effective radius, and the position of the centre. We have two 
sets of results, depending on how the fit was performed (see 
(iv) below). The software we developed works as follows: 

(i) The code generates a 2D projected image of the galaxy 
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Figure 4. Comparison of photometric parameters from the 
isophotal profile fit (x-axis) with those from one of the 2D fits 
(the least x^, y-axis). Subscripts 6 and 8 indicate filters F606W 
and F814W. 



using the r^'''* law on a subsampled grid. The pixel size 
is chosen to be approximately ten times smaller than the 
galaxy effective radius in order to make discretisation prob- 
lems negligible. 

(ii) The resulting "ideal" model is rebinned to the 
WFPC2 chip scale. 

(iii) The rebinned model is convolved with the Tiny Tim 



( |Krist 1994| ) PSF. The Tiny Tim (non resampled) PSF in- 
cludes the effects of diffraction (i.e. of the telescope and cam- 
era optics) and the Pixel Response Function (i.e. the spread 
caused by electron diffusion in the WFPC2 CCDs) . 

(iv) This model is fitted to the data either by least 
or by least squares, using a simplex algorithm to find the 
minimum (amoeba, Press et al. 1992). 

For each galaxy we ran the code twice, minimising both 
the and the unweighted residuals (least squares). In the 
case of G the fit was not straightforward because of the 
companion H. Therefore we subtracted from the data, be- 
fore performing the fit, a model for H, obtained from an 
r^/"* profile with re, SBe, pa, and e taken from the isophotal 
profile fit. The result was stable with respect to variations of 
the parameters of H. In Figure I and I the original WFPC2 
image and the fit residual are shown for all the galaxies fitted 
with the 2D fits. 



3.2.3 Results 

All galaxies are well fitted by the r^^* law (de Vaucouleurs 
1948) with the exception of galaxy B which shows face-on 
spiral-like residuals and might be an early-type spiral. 

Standard magnitudes ar e obt ained from the data fol- 
lowing the Holtzman et al. ( 1995 ) calibration. In Figure ^ 
we plot the photometric parameters obtained from one of 



the 2D fits (the least x ) versus those given by the isophotal 
profile fit. The results are in very good mutual agreement, 



Table 3. Photometric parameters through filter F606W. Errors 
(see Section 3.4) do not include sky-subtraction contributions, 
which we estimate to be 0.5 per cent on re, 0.05 mag arcsec"^ 
on SBe and 0.02 mag on m (see text for details). 



galaxy 


re 




ooe 








A 


1.65 


0.13 


20.64 


0.13 


17.57 


0.12 


B 


0.588 


0.038 


19.73 


0.11 


18.90 


0.05 


D 


0.793 


0.042 


20.99 


0.09 


19.5 


0.03 


E 


0.794 


0.071 


21.60 


0.13 


20.12 


0.06 


F 


1.145 


0.017 


22.56 


< 0.01 


20.27 


0.03 


G 


1.589 


0.064 


21.81 


0.07 


18.81 


0.01 


I 


0.636 


0.046 


20.92 


0.13 


19.92 


0.04 


O 


0.795 


0.011 


23.15 


0.05 


21.65 


0.02 


P 


0.526 


0.029 


22.20 


0.11 


21.66 


0.01 


Q 


0.612 


0.018 


23.20 


0.06 


22.27 


0.02 



Table 4. Photometric parameters through filter F814W (same as 
Table 2). 



galaxy 


rc 


Src 


SBe 


(5SBe 


m 


(5m 


A 


1.99 


0.15 


19.92 


0.19 


16.44 


0.13 


B 


0.548 


0.014 


18.58 


0.05 


17.89 


0.04 


D 


0.645 


0.058 


19.29 


0.16 


18.27 


0.04 


E 


0.872 


0.080 


20.31 


0.12 


18.64 


0.07 


F 


1.29 


0.12 


21.33 


0.11 


18.81 


0.10 


G 


1.419 


0.075 


20.29 


0.09 


17.53 


0.02 


I 


0.606 


0.051 


19.47 


0.15 


18.58 


0.03 


O 


0.654 


0.007 


21.02 


0.05 


19.94 


0.02 


P 


0.502 


0.018 


20.48 


0.08 


19.98 


< 0.01 


Q 


0.547 


0.053 


21.14 


0.17 


20.48 


0.03 



even though the procedures differ in some important details, 
e.g., the isophotal profile fit is in magnitudes and allows for 
pa and e variation, while the 2D fits are in counts and have 
fixed pa and e. As final value of the parameters we use the 
average of the three results; the photometric parameters are 
shown in Tables ^ and ^. Integrated luminosities are those 
of the best-fitting r^^* law. 



3.3 A cross-check with the growth curve 
technique 

The growth curve technique is a standard procedure for 
nearby samples. Two problems of this technique are: i) 
the statistical dependence between different aperture mag- 
nitudes, e.g. in the presence of light coming from a nearby 
galaxy or of diffuse background sources, the systematic error 
affects the flux within all the apertures; ii) the sensitivity 
to sky-subtraction errors. We measured the photometric pa- 



rameters by fitting the growth curve of an 



1/4 



profile to 



the aperture magnitudes and the results were in reasonably 
good agreement with the other two methods. 



3.4 Error analysis 

We define as internal error of the parameters the standard 
deviation of the results of the two 2D fits and the isopho- 
tal profile fit. Since photon shot noise is negligible in our 
context, the main source of error not included in the inter- 
nal error is the uncertainty due to sky subtraction (0.5 per 
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Figure 2. WFPC2 images and 2D fit residuals for the galaxies with measured kinematics: A,D, E, F, G, I. Prom the left to the right 

wc display the F606W image, the F606W residual, the F814W image and the F814W residual. The images (and their sizes) arc the ones 
used for the 2D fits, see Table 1. Galaxies A, D show a residual in the innermost pixels which is in agreement with the deviation from 
the r^/** law in the luminosity profile. The bad column in the centre of E is clearly evident in the residual. F shows a peculiar pot handle 
structure just to the upper left of the bad column. 
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Figure 3. WFPC2 images and 2D fit residuals for the galaxies without measured kinematics: B,0,P,Q. In the residual of B a spiral 
pattern is recognizable. 



cent on re, 0.02 mag on m, 0.05 mag arcsec"^ on SBe). We 
estimated this uncertainty by measuring how the photomet- 
ric parameters changed with the sky level, and we added it 
quadratically to the internal error, so as to obtain the total 
error; more precisely we used the variation of the photomet- 
ric parameters when the sky is shifted by one standard de- 
viation, calculated as the scatter of the sky levels measured 
in different areas (a typical value is 2%). The photomet- 
ric parameters re and (7) e are correlated, as can be seen 
from the errors calculated directly and reported in Table ^ 
The particular combinations of effective radius and surface 
brightness are chosen because they enter the definition of 
the K space (Bender, Burstein & Faber f992). Also the er- 
ror on the colour (mpgQg-^y — 'rn■pgl4^J\l) is different from the 
sum of errors on the single filter magnitudes (see Table ^ 
since the errors in the two filters are correlated: it generally 
happens that when one method (e.g. 2D x^) gives a lower 
magnitude than another (e.g. profile fits), the same thing is 
repeated in both filters. 



Table 5. Photometric parameters II. Errors are shown as ex- 
amples of the correlation of the photometric uncertainties. The 
quantity colour is defined as 'TipgOGW ~ '"F814W' 



jalaxy 


colour 


(5colour 




<5(/>ere/(/)er-e 


A 


1.127 


0.012 


0.23 


0.089 


B 


1.015 


0.012 


0.24 


0.047 


D 


1.226 


0.018 


0.19 


0.027 


E 


1.486 


0.013 


0.29 


0.034 


F 


1.463 


0.068 


0.025 


0.011 


G 


1.273 


0.007 


0.15 


0.028 


I 


1.333 


0.005 


0.28 


0.044 


O 


1.705 


0.006 


0.091 


0.028 


P 


1.678 


0.013 


0.26 


0.045 


Q 


1.789 


0.051 


0.15 


0.031 



4 SPECTROSCOPY 

The spectroscopic observing run took place in the period 
April 18-22 1996 at the 3.6-m ESQ telescope. Long slit (I'.'S 
slit width) spectroscopy was obtained with the EFOSC spec- 
trograph, with the Orange 150 grism in the spectral range 
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Table 6. Spectroscopic data: number of observations (nexp), to- 
tal exposures times (texp), and average signal-to- noise ratio per 
pixel (S/N). 



galaxy 


ncxp 


texp (s) 


S/N 


A 


4 


vonn 
( zUU 


DD 


B 


4 


7200 


40 


D 


4 


12000 


22 


E 


5 


13200 


17 


F 


5 


13200 


12 


G 


3 


10800 


35 


I 


3 


10800 


22 


LM 


2 


4800 




N 


2 


4800 




O 


10 


36000 


7 


P 


10 


36000 


6 


Q 


10 


36000 





between approximately 5200 and 7000 A (512x512 CCD, 
with pixel size 30 equivalent to about 0'.'57). The slit 
length of about 180" projects to only about 320 pixels on 
the CCD. For each target we obtained multiple exposures. 
The estimated seeing, measured with R-band acquisition im- 
ages, was O'.'S FWHM. Table ^ gives the basic observation 
log. 

The slit positions were carefully chosen to include as 
many objects (two or three at a time) as possible. Every 
night we obtained an exposure with the He-Ar lamp at 
zenith; before and after every galaxy spectrum we took a He 
lamp spectrum as a check of the dispersion relation shifts. 
We observed two spectrophotometric standard white dwarfs 
in order to perform relative-flux calibration of the spectra. 
The identification of our targets as normal, non star-forming, 
ellipticals was confirmed by the absence of emission lines in 
their spectra. 

4.1 Reduction 

A superbias file was obtained by averaging about a hundred 
bias images and was subtracted from the data. Dome flat- 
field exposures were used to produce a flatfield along the 
wavelength direction (CCD columns). The slit function was 
obtained from sky frames and was found to be independent 
of wavelength. Dark current subtraction was done by pro- 
ducing a dark frame from the average of 8 darks 1 hour 
long. Since this mean dark did not contain structures, we 
only subtracted the mean value. 

Throughout the reduction and analysis we have avoided 
(non-linear) rebinnings. To this end, we did the reduction 
before wavelength calibration. The extraction of the spec- 
tra was a crucial point of the reduction process, given the 
faintness of the signals. Averaging too many rows would have 
meant too much noise, using only the central row would have 
meant losing important signal. As we could not simply ex- 
tract a fixed pixel width (along the spatial direction) because 
of optical distorsions, we used a code developed by M0ller 
& Kjfflrgaard (1992) to obtain a one-dimensional spectrum, 
extracted with optimal signal-to-noise weighting. The code 
also removes cosmic rays, identified as pixels five standard 
deviations away from the expected signal. All spectra were 
carefully inspected by eye before and after extraction to 
check the proper cosmic ray removal. 



Given the faint signal of the spectra, we took special 
care during sky subtraction. For each galaxy we extracted 
two sets of columns, to the left and to the right of the galaxy. 
The columns were selected to be as close as possible to the 
galaxy to limit optical distorsion problems. Columns of one 
set were combined neglecting the highest and lowest pixels 
to remove cosmic rays. The result was eye-inspected and 
then the two spectra obtained were averaged to obtain a 
high signal-to-noise sky spectrum. This procedure allowed 
us to check the optical distorsion by comparing the emis- 
sion lines between the left and the right spectra. The sky 
spectra obtained had broader lines than the single column 
ones, because we combined spectra wavelength-shifted by 
optical distorsion. Sky contribution was subtracted during 
spectrum extraction by deconvolving the artificial additional 
width. The resulting one-dimensional spectra were combined 
to obtain the final galactic spectrum, after verifying that the 
relative shifts of the positions along the y-axis of the emis- 
sion lines were negligible. The reduced spectra are shown in 
Figure ^. The average signal-to-noise ratio per pixel of the 
spectra is listed in Table The average signal-to-noise ra- 
tio ranges from 65 per pixel for the brightest galaxy (A) to 
12 per pixel for the faintest galaxy for which we obtained 
a velocity dispersion (F); the estimated error is accordingly 
larger for F than for A. For comparison, Davies et al. (1987) 
in their low redshift study, with a total signal-to-noise ra- 
tio equivalent to 5T0'' photons, estimate the random error 
to be about 10 %. The analysis in Section ^ is done both 
with the entire sample and excluding the two galaxies E and 
F with the lowest signal-to-noise spectra among those with 
measured velocity dispersion. 

4.2 Wavelength calibration and measurement of 
the instrumental resolution 

We found the dispersion relation for each He-Ar lamp us- 
ing the MIDAS GUi/long software. The different calibrations 
were all within the errors. We used short He lamp images to 
check that shifts between single exposures were negligible. 
The He-Ar lamps were also used to measure the instrumental 
resolution: the line shapes were fitted with Gaussian func- 
tions to obtain line widths as a function of CCD pixel num- 
ber, and, via dispersion relation, of wavelength. The FWHM 
is independent of A within the measurement scatter and it 
was determined to be: 

AA = (8.39±0.26)A, (2) 

Translated into velocities it corresponds to a kinematic res- 
olution ranging from 206 to 155 km s~^ from the blue to the 
red. The 3 % error is due to the scatter in the measurement 
(see also Section 4.4) . The instrumental resolution measure- 
ment is very important for the kinematic fit, as we discuss 
in Section 4.3. Therefore we carefully examined all lines, to 
avoid ghost lines which were found and which would have af- 
fected the measurement. We checked the resolution by mea- 
suring the width of the sky lines and the results were con- 
sistent with the ones found with the He-Ar lamp (8.11±0.39 

A). 

For comparison, the instrumental resolution of the 
Lick/IDS spectrograph was between 8-10 A at shorter wave- 
lengths (see e.g. Dalle Ore et al., 1991). The velocity disper- 
sions we find are larger than half the instrumental resolution. 
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which is commonly considered as a limit (see e.g. Dressier 
1979; Kormendy 1982; Bender, 1990). It has been noticed 
that when the velocity dispersion becomes of the order of 
the kinematic resolution and the signal-to-noise is low a sys- 
tematic error may be introduced (see e.g. j0rgensen, Franx 
& Kjffirgaard, 1995). In order to study this effect we broad- 
ened some stellar templates to 150, 200, 250, 300 km s~^and 
we added noise, thus producing artificial galaxies with S/N 
per pixel ranging between 10 and 60. The artificial galaxies 
were then fitted with the same stellar templates to mea- 
sure how well the original velocity dispersions were recov- 
ered. The largest velocity dispersions were recovered within 
few percents (2-4 %) down to S/N=10, while the smaller 
one (150 km s"'^) resulted systematically underestimated 
by ~ 10% at S/N 10. No systematic trends were noticed 
above S/N=15. This bias is negligible to our measurement, 
as galaxies E and I, the only with measured velocity dis- 
persion below 200 km s~^(see Table 7), have high enough 
S/N (see Table |), while F, the only with S/N below 15 has 
a ^ 200 km s"^ 



4.3 Kinematics 

In order to derive a robust measurement of the kinematics 
of our sample of galaxies we used two different and inde- 
pendent cod es: i) the Gauss-Hermite Fourier Fitting 



Software (Franx, Illingworth fc Heckman 1989| ; 



van del 



Marel fc Franx 1993), hereafter GHFF, and ti) the Fourier 



Quotient (Sargent et al. 1977), hereafter FQ, in the version 
modified by Rose (Dressier, 1979) and used e.g. by Stiavelli, 
M0ller & Zeilinger (1993). The results of the two methods 
are mutually consistent, within the estimated errors. How- 
ever, the GHFF, allows a reliable estimate of the errors (see, 
e.g., Rix & White, 1992, and our discussion below) while 
the version of the FQ we used does not; moreover methods 
like the GHFF are known to be less sensitive to template 
mismatches than the FQ (see e.g. Bender 1990, and our 
discussion below). Furthermore the GHFF provides more 
insight in the fitting procedure, because the continuum sub- 
traction, the tapering, and the filtering of both observed 
and template spectrum can be checked step by step. Finally 
it provides a quality parameter (the x^) Sind the residual 
spectrum. Therefore we shall consider the values obtained 
with the FQ only as a check of the values obtained with the 
GHFF. We report the complete results obtained with both 
codes (Tables 7 and 8 for the GHFF, Table 9 for the FQ). 
The discussion of the FP properties in Sections ^ and ^ is 
based on the GHFF method, but the FQ method leads to 
the same conclusions. 



4-3.1 Template spectrum 

We used star spectra from the Jacoby et al. ( 198^ ) library as 
starting bricks for the template construction. We simulated 
the measured experimental resolution, taking into account 
the non-negligible width of the original Jacoby spectra (4.5 
A FWHM), by convolving the redshifted template spectra 



z is the galaxy redshiftQ. The choice of the correct spectral 
class for the template is critical in order to avoid system- 
atic errors. Therefore, we decided to fit the galactic spectra 
with each of the spectral types from FO to K7 to be sure 
to include all the most relevant spectral types. By setting 
some fit quality criteria (see the following section) we se- 
lected a set of "good" spectral types. The range of variation 
of the central velocity dispersion within the "good" spectral 
types provides a check of the robustness of the result, and at 
the same time an estimate of the systematic error possibly 
induced by residual template mismatches. We avoided the 
alternate approach of using synthetic galaxy spectra because 
of its heavy model dependence on the underlying stellar pop- 
ulation. 



4.3.2 The Gauss-Hermite Fourier Fitting 
Software ('GZ/FF; fit 

The GHFF, developed by R. P. van der Marel & M. Franx 
(Franx, Illingworth & Heckman, 1989; van der Marel & 
Franx, 1993) models the broadening function B as the first 
terms of a Gauss-Hermite series, thus allowing one to derive 
a better estimate of the velocity dispersion in case of non- 
Gaussian broadening. The advantages and caveats of this 
approach are discussed in detail in van der Marel & Franx 
(1993). In order to achieve a meaningful estimate of the first 
non- Gaussian terms however, very high signal-to-noise ra- 



Bender 



with a gaussian of width at = 8.39^ - (4.5(l+z))^ where 



tio is requir ed (Franx, Illingworth & Heckman 198!; 
et al. 1994) and generally the studies on the Fundamental 
Plane with this software fit only the Gaussian component 
(see e.g. j0rgensen, Franx & Kjaergaard, 1995). Therefore 
we used the code in Gaussian- fitting mode. 

The GHFF code provides the redshift (z^gj^g-), the line 
strength ratio (7giiff) and the velocity dispersion (cgiiff). 
The fit gives also an error value for the parameters that we 
will call Sv^Yiff' hghS ' '^'^ghff- 

Generally, the templates within the spectral types GO- 
K2 provided the best fit, as could be judged by looking at the 
X^, the residuals and the broadened template, but no tem- 
plate was clearly a better approximation than the rest. In 
order to better understand the systematic trends, we created 
some "artificial galaxies" by broadening the templates and 
adding noise to reproduce a signal-to-noise ratio of about 
12 per pixel. The GHFF with the templates in the range 
G0-K2 was able to recover the correct velocity dispersion, 
while generally the F stars gave lower sigmas and the K3-7 
higher ones. The scatter in the values obtained with differ- 
ent noise patterns was comparable with the formal error of 
the fit when the same star was used both as "galaxy" and 
template, but was higher even in case of small mismatch- 
ing. For example a "galaxy" produced and fitted with a G5 
spectrum gave on average 96 % the correct value, with 15 % 
average formal error and a measured scatter of 17 %; using a 
G2 star as template we obtained equally 96 % of the correct 
value, with 22 % scatter and 17 % average formal error; with 



t The redshift was first determined by recognizing the most 
evident absorption lines and then by iteratively refining the 
measurement. 
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Figure 5. Wavelength (in A) calibrated spectra of the galaxies with measured velocity dispersion. The shaded spectral regions were 
aflfected by strong sky-emission lines. These regions have not been used in the kinematic fit. The standard deviation per pixel is also 
shown (lower curves). The arrows identify some of the main absorption features, such as Ca H (3968.5 A) and K (3933.7 A), the G band 
(4304.4 A), Mgi, (5175.4 A) and Na D (5892.5 A). 



a G6 the same figures were 104 %, 21 % and 19 %. There- 
fore we decided to keep only the G0-K2 interval, where no 
systematic trend was noticed. 

We defined the scatter of the velocity dispersion of the 
"good" templates to be the estimate of the systematic er- 
ror due to residual mismatches. In order to obtain results 
comparable to the FQ ones, we selected the good fits with 
criteria applicable also to the FQ (that in the version we have 
does not provide the x^)- We used the information provided 



by the code as an indicator of the quality of the fit. Specif- 
ically, we used the line strength ratio Tg^ff as an estimator 

of the mean ratio of the abundances of the absorbers (i.e., 
roughly speaking, the metallicity) and Sag^ffas a measure 
of the mismatch in the shape (primarily due to differences 
in the spectral type) . A fit is considered acceptable if it sat- 
isfies the three conditions: i) 0.5 < 7giiff < 1-5, ii) Ugj^ > 
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Table 7. Kinematic results from the GHFF fit I. For each galaxy 
we list the average central velocity dispersion (o"g]-|ff ) together 
with the random error {S ""gjiff ) ^.nd the systematic error (A 
(jgjjg ) • "^^^ random error is the fit formal uncertainty for the 
best-fitting template, while A o"gjjff is the scatter of the values 
obtained with the "good" templates (see text). 



gal 


'^ghff 


^o-ghff 


^'^ghff 


A 


232 


11 


28 


D 


271 


24 


16 


E 


175 


28 


28 


F 


195 


41 


74 


G 


226 


15 


25 


I 


166 


21 


29 



Table 9. Kinematic results from the FQ fit. For each galaxy 
we list the average central velocity dispersion, redshift and line 
strength ratio ((T£q ,2£q , 7£q ) together with the systematic error 
(A (T£q , A Z£q , A 7£q ), that is the scatter of the values obtained 
with the "good" templates (see text). 



gal 


'^fq 


^'^fq 


^fq 


^^fq 


Tfq 


A7fq 


A 


208 


39 


0.1462 


MO-4 


0.92 


0.16 


D 


252 


13 


0.3844 


1-10-4 


0.84 


0.14 


E 


211 


10 


0.2935 


2-10-4 


0.67 


0.10 


F 


233 


74 


0.2934 


2-10-4 


0.61 


0.01 


G 


216 


55 


0.2924 


1-10-4 


0.56 


0.03 


I 


214 


39 


0.2918 


2-10-4 


0.65 


0.10 



Table 8. Kinematic results from the GHFF fit II. Same as Table 
^ for the redshifts {Zg\^ff) and line strength ratios (7gjifi^). 



gal 


^ghff 


<5 ^ghfl^ 


A 


^ghff 


Tghff 


^ Tghff 


^ Tghff 


A 


0.1462 


310-5 


3 


10-4 


0.90 


0.02 


0.14 


D 


0.3844 


5-10-5 


1 


10-4 


0.80 


0.03 


0.04 


E 


0.2936 


6-10-5 


1 


10-4 


0.91 


0.04 


0.14 


F 


0.2933 


1-10-4 


1 


10-4 


0.95 


0.10 


0.19 


G 


0.2942 


4-10-5 


1 


10-4 


0.93 


0.03 


0.16 


I 


0.2924 


4-10-5 


1 


10-4 


0.84 


0.04 


0.16 



and til) (^fghfj <o"giiff . As final value we adopt the average 
of all acceptable results weighted on the fit formal errors. 

We checked what templates were rejected by this "x^ 
blind" algorithm, and the criteria were very effective in se- 
lecting the best a-nd the best matching templates. The 
GHFF compared to the FQ produced generally a larger 
number of fits within the quality criteria and a lower scatter 
as a function of template, thus confirming to be less sensitive 
to template mismatches. As a further check we also consid- 
ered a simple average of the five best a-nd the changes 
were much smaller than the error bars. The average values 
are listed in Table ^ The fit formal errors listed refer to the 
template which gave the least x^- 

We have omitted galaxy B in our kinematic discussion 
since its kinematic fit does not satisfy our quality criteria 
and indicates a very small velocity dispersion, below our 
resolution limit and consistent with the morphological ap- 
pearance as a face-on early spiral (see Figure ^). 



4.3.3 Fourier Quotient fFQj /i< 

The basic idea of the FQ is the hypothesis that one can de- 
scribe a galaxy spectrum G as the convolution G = B o S 
of an average star template spectrum S with a broadening 
function B. The broadening function is purely kinematic, 
while the template spectrum includes the broadening due to 
the instrumental resolution. The template spectrum has to 
be a good approximation of the galaxy spectrum especially 
in its absorption features because FQ works on the spectra 
after subtraction of the continuum. Selecting the best tem- 
plate is cri tical t o the measurement since, as shown, e.g., 
by Bender (1990), FQ is very sensitive to the spectral type 
of the star used as a template. To perform the kinematic fit 
with the FQ we needed spectra free of their continuum com- 



ponent. The continuum was subtracted out by interpolating 
between regions with no spectral features. 

The version of the FQ code available to us requires a 
template and a galaxy spectrum linearly rebinned in wave- 
length and then it rebins them in log A. We modified it to 
avoid the extra rebinning. The new version requires the one- 
dimensional galaxy spectrum as a function of the CCD pixel 
number and the dispersion relation given as polynomial co- 
efficients; the rebinning is done directly in log A. The spec- 
tral regions affected by the strongest sky emission lines were 
masked out. 

Free parameters in the fit are the redshift (zfq), 
the width (fj-q), and a normalising factor (7fq , the line 
strength). The fit gives also a formal error value for the pa- 
rameters that we will call 5v^^, (57fq , 5o"fq. 

A detailed study of the results shows no correlations 
between fit parameters and errors. The only general trend 
is that best results are obtained for templates of spectral 
type G, and they become worse towards the edges of the 
range we considered, FO and K7 (for the complete results 
see Treu 1997), consistently with the findings of the GHFF. 
Generally, the F star templates gave systematically lower 
velocity dispersion, while the K3-7 templates gave system- 
atically higher values (see also the discussion in the GHFF 
section). In Table ^ we list the average velocity dispersion, 
line strength ratio and redshift, together with the scatter in 
the results between different templates as for the GHFF. 

Consistently with what is found with the GHFF, the 
velocity dispersion derived for B was below our resolution 
limit. 



4.4 Error analysis 

The errors on a, 7, z are a combination of a systematic com- 
ponent due to template errors and mismatch, and a random 
component due to the limited signal. For the systematic er- 
ror we adopt the standard deviation Act of the sample of 
the acceptable a values (see Tables ^ and ^). For the ran- 
dom error, we adopt the formal error given by the GHFF 
fit with the best-fitting template, although it is probably an 
overestimate because it includes possible residual template 
mismatches. The random and systematic errors for the line 
strength ratio (7) and the redshift (z), defined in the same 
way, are hsted in Tables § (FQ) and | (GHFF). 

The uncertainty in the measurement of the instrumental 
resolution introduces a systematic error in the measurement 
of the velocity dispersion that is in our case negligible. For 
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example, considering A, the galaxy with the highest signal- 
to-noise spectrum, the random error is estimated by us to be 
5 %, the systematics due to template mismatches 12 %, and 
therefore a further 3 %, to be added in quadrature, does not 
change the total error significantly. As a double check we re- 
peated the fits with templates with different resolution {Sat 
equivalent to 0.26 A) and the results changed as expected: 



Table 10. Galactic extinction coefficients. 



(3) 



For galaxy I we obtained the larger difference (7 km s~^), 
still negligible with respect to the other sources of error. 

The spectrum of galaxy D includes the lines Ca H and 
K. While widely used to perform such measurements (e.g. 
Dressier 1979), these lines has been reported to induce a 
slight overestimate of the velocity dispersion (Kormendy, 
1982, estimates it to be between 5 and 20 %; see also Ko- 
rmendy and lUingworth, 1982). They suggest the problem 
may be caused by the intrinsic width of the lines or by the 
steepness of the continuum in that spectral region. Ifowever, 
Dressier (1984), suggest that they may be best suited for 
measurement of large velocity dispersion for faint objects, 
such as D. For these reasons we performed the kinematic fit 
also exluding the region of Ca H and K and we found a ve- 
locity dispersion 10 % higher in this case, a variation of the 
order of the estimated errors. In the analysis presented in 
this paper we shall consider the results from the entire spec- 
trum, for its higher signal-to-noise ratio. Nevertheless, when 
a larger sample of these objects will be available we plan to 
better address this issue by studying systematic variations 
induced by Ca H and K on the sample. 

Similarly the spectrum of A includes the region of NaD, 
which may be affected by interstellar absorption even in ell- 
pitical galaxies (e.g. Dressier, 1984). We repeated the anal- 
ysis exluding the NaD region and found that this resulted 
in a 3% increase in the velocity dispersion, again within the 
estimated errors. 

The error on z is caused by uncertainties in the wave- 
length calibration (the dominant source), in the fit results 
and in the absolute motion of the template stars, for the 
Jacoby et al. (1984) spectra were not corrected for peculiar 



motions (estimated to be of order 0.0001 by the authors). 
The total error adds up to about 0.0005. In the cases with 
insufficient signal, the redshift has been measured by iden- 
tifying the main spectral features, with a total estimated 
error of ~ 0.005. These redshifts were then been checked 
via superimposition of redshifted templates. The listed er- 
ror corresponds to ~ 10 pixel in the studied spectral range. 

4.5 Relative flux calibration 

The instrumental throughput is not constant as a function 
of wavelength. In order to obtain a good relative flux cali- 
bration for our spectra, we observed two spectrophotometric 
standard white dwarfs through our instrumental setup for a 
total of five exposures. By comparison with the spectra taken 
from the MIDAS libraries, we obtained the response function 
separately for each exposure, using MIDAS GUi/long a nd th e 
values of galactic extinction from Burstein & Heiles ( 1982 ). 
The results were all within two percent and we adopted their 
average function to correct the spectra for the instrumental 
sensitivity. 



field 



E(B-V) ApgQg-\Y ApsUW 



u5405 


0.18-0.21 


0.57 


0.37 


ur610 


0.03-0.06 


0.13 


0.09 


urzOO 


0.06-0.09 


0.22 


0.14 


ustOO 











utSOO 


0.03-0.06 


0.13 


0.09 



5 PHOTOMETRIC AND KINEMATIC 
CORRECTIONS 

The spectroscopic and photometric results described so 
far require a number of corrections. Photometry should 
be corrected for galactic extinction (Subsection 5.1), while 
WFPC2 magnitudes should be t ran sformed to standard 
restframe magnitudes (Subsection 5.2) and fg^^g- should be 
correct ed in to a standard central velocity dispersion a (Sub- 
section 5.3). 



5.1 Galactic extinction 

For galactic extin ction , we used E(B-V) values from the 
Burstein & Heiles ( |l982| ) maps, and the relations ^F606W ~ 
2.907 E(B-V ) and Ap^i^yj^ = 1.902 E(B-V), calculated by 
Romaniello (1997). The extinction coefficients obtained for 
the various fields are given in Table hol 



5.2 An improved way to calculate K-correction 

The photometric parameters measured through WFPC2 fil- 
ters have different physical meaning depending on the galaxy 
redshift. The common technique to obtain standard magni- 
tudes is to introduce the so-called K-correction, 



K; 



-2.5 log 



/ Fx{\)So{X)dX 



J Fx{X)So(X{l + z))dX 



(4) 



where Fx is the spectral flux density of the galaxy and So 
is the response function of the instrumental setup used dur- 
ing the observations. To calculate the K-correction at in- 
termediate redshift, knowledge of a large spectral range is 
required. In order to reduce the required spectral range and 
to minimise the correction, it is convenient to calculate the 
restframe absolute magnitudes (Afw) in different filters from 
the observational ones (mo): 



Mw = Amwo{z; Fx) + mo - DM{z;Q.,Q.a), 
where 



(5) 



Amwo 



-2.5 log 



/Ff (A)Sw(A)dA 
/ F^(\)So{Ml+z))d\ 



+2.5 log 

DM{z;Q,nA) is the bolometric distance modulus. 



(6) 



Fx^ is the spectral flux density of a standard star 



and 
For 

each redshift, we calculated corrections from F606W to 
Johnson B (AmgpgQg-^^) and from F814W to Johnson V 

(AmvF814w)- 

To perform this calculation we chose the following, 

purely empirical, approach. For each galaxy we selected from 

the synthetic library of Bruzual & Chariot (1993; GISSEL96 
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Table 11. Photometric corrections 
galaxy z 



Am 
B F606W 



Am 
V F814W 



A 


0.146 


1.11 


1.12 


B 


0.146 


1.01 


1.04 


D 


0.384 


0.44 


0.86 


E 


0.294 


0.68 


1.00 


F 


0.293 


0.68 


0.99 


G 


0.294 


0.63 


0.92 


I 


0.292 


0.65 


0.94 


O 


0.647 


-0.28 


0.65 


P 


0.550 


-0.05 


0.74 



For the galaxies B, O, P with insufficent signal-to-noise ratio, the 
redshift was found by identifying the strongest spectral features. 
The redshift of P is somewhat less certain because the signal-to- 
noise ratio is particularly low. 

version) the spectrum which best approximates the galac- 
tic spectrum under consideration in the observed spectral 
range. Assuming that the selected spectrum provides a good 
description of the spectral energy distribution of our galaxy 
also outside the observed wavelength range, we used it to 
compute Amwo- This method also provided us with an er- 
ror estimate based on the range Amwo obtained from all 
the spectra not significantly different from the observed one; 
the typical error on the K-correction is ~ 0.05 mag. Since we 
were particularly interested in the continuum shape, which 
is the main contribution to large-passband fluxes, we needed 
flux c alibrated spectr a. The so-called age-metallicity degen- 
eracy ( Worthey 1994) with respect to the continuum shape is 
not a significant concern, because we were only choosing the 
best phenomenological approximation. In Table |l^ we list 
the computed Amwo values. As a cross-check we compared 
the WFPC2 colours of our galaxies with the WFPC2 colours 
"measured" on the redshifted synthetic spectra. The colours 
from the best approximation spectrum were all within 0.05 
magnitudes from the observed one. 



5.3 Reduction to a standard central velocity 
dispersion 

In Section I.3.3| o-gj^g-was defined as the broadening of the 
projected spectrum integrated over the solid angle sub- 
tended by the CCD. The velocity dispersion in elliptical 
galaxies varies with distance from the center. Usually, the 
velocity dispersion is defined as the kinematic broadening 
of a spectrum integrated over a reference solid angle. The 
choice of a standard reference solid angle makes it possi- 
ble to compare values obtained in different studies. Previ- 
ously used reference solid angles are, e. g., a fixed angular 
size at a fixed distance or a fraction of the effective radius 
(J0rgensen, Franx & Kjargaard 1996; hereafter JFK96). 
The latter choice is to be preferred because it is related to an 
intrinsic length scale of the galaxy and does not depend on 
cosmological parameters. Following JFK96, we correct all 
measurements to a circular aperture of re/4. To calculate 
the correction we need to model the large-scale kinematic 
profile of elliptical galaxies. A power law: 



aMoc(^)', 



Table 12. Kinematic correcting factors B with standard devi- 
ations SB (see text) and final central velocity dispersion cr in 
km s~^ . 



galaxy 



'^ghff 



B 



SB 



A 


232 


1.07 


0.04 


248 


D 


271 


1.14 


0.08 


309 


E 


175 


1.10 


0.06 


192 


F 


195 


1.08 


0.05 


211 


G 


226 


1.09 


0.05 


246 


I 


166 


1.11 


0.06 


184 



is generally a good description of the kinematic profile (see, 
e.g., the radially extended profiles in CaroUo & Danziger 
1994a, 1994b, Bertin et al. 1994 with -0.1 < d < 0). 
We assume that 



2-Krdra {r)DV{r) 



(8) 



where r is the projected angle variable, A is the projection 
on the focal plane of the solid angle subtended by the instru- 
ment and DV{r) is the r^''* law, appropriately normalised. 
Under these assumptions 



2 lo 
^ '^Shff— 



1/4 



dxxa^{xre)DV{xre) 
J^drra'2{r)DV{r) 



(9) 



(7) 



The correcting factors were computed by taking into account 
the number of pixels used during each spectrum extraction 
and are summarised in Table Given the lack of informa- 
tion about the kinematic profiles, the correcting factors is 
taken to he B = [S(— 0.1) -I- B{0)]/2, with an error estimate 
SB = [B{-0.1) - 5(0)1/2^3 (-5%). 



5.4 Effective radii 

The radii RqY ^^'^ -^eB calculated for the cen- 

tral rest wavelength of each filter by linear interpola- 
tion/extrapolation in wavelength between the measured 
radii in the WFPC2 filters F606W and F814W at their ob- 
served central wavelengths. We take the central wavelengths 
of B, V, F606W and F814W to be 4400, 5500, 5935, and 
7921 A, respectively, resulting in the following explicit for- 
mulae for the angular sizes: 

res = {r-eFfi()fiw[7921 - (1 + 2)4400] 

+reF8i4w[(l + 2)4400 - 5935]}/1986 (10) 
rev = {reF606w[7921 - (1 + 2)5500] 

+reF8i4w[(l + 2)5500 - 5935]} /1986 . (11) 
The typical correction is of order of 3%. 

6 RESULTS 

We now consider the astrophysical implications of our mea- 
surements in terms of the changes in scaling laws as a func- 
tion of redshift. We emphasize that our sample, although 
small (only six galaxies with a reliable measurement of the 
velocity dispersion), is unique in that it focuses on field 
galaxies. Therefore, we will first consider our sample by itself 
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Table 13. Final corrected photometric values. Effective radii are 
expressed in h^Q kpc. B and V refer to the corresponding Johnson 
passbands. Quantities are computed using Q = 1, Q/^ = 0. 



galaxy 




SBeB 




^eV 


SBeV 


Afy 


A 


5.02 


21.03 


-21.23 


5.73 


20.37 


-22.30 


B 


2.02 


20.02 


-20.00 


1.94 


18.95 


-20.94 


D 


4.93 


19.89 


-22.16 


4.22 


18.65 


-22.94 


E 


4.30 


20.60 


-21.13 


4.60 


19.83 


-22.09 


F 


6.18 


21.57 


-20.97 


6.76 


20.85 


-21.92 


G 


9.02 


20.75 


-22.5 


8.39 


19.72 


-23.29 


I 


3.52 


19.89 


-21.35 


3.40 


18.93 


-22.20 


O 


5.61 


20.70 


-21.83 


4.58 


19.50 


-22.61 


P 


3.58 


20.29 


-21.21 


3.43 


19.32 


-22.09 



(Subsection 3.1), and then proceed to compare with clust er 
galaxy measurements from the literature (Subsection 6.2). 



6.1 Results from our sample: Scaling laws for field 
galaxies at intermediate redshifts 

The input data for the galaxies in our sample are presented 
in Table for the central velocity dispersion and in Table 
for the photometry. Absolute magnitudes and metric sizes 
in Table ^ are computed using Hq—SO hso km s^^Mpc"^, 
n = 1, = 0. 

We consider tw o scaling laws: t he FP (Eq. ^ and the 
Kormendy relation (Kormendy 1977) 



SBe = Hi log Re +a2. 



(12) 



The latter is useful for our purposes because it does not 
involve internal kinematics, and thus it can be followed to 
higher redshifts {z = 0.647 for our sample), although it has 
a larger intrinsic scatter than the FP and its coefficients 
are much more sensitive to sample selection biases (see e.g., 
Capaccioh, Caon & D'Onofrio 1992). 

Because of the small number of points in our sample, 
we do not try to rederive the parameters a and /3 for the FP 
independently, but we compare our sample of field ellipticals 
with the local relations, adopting their slopes. In Figures ^ 
and ^ we compare our galaxies with a sample of Coma el- 
lipticals in Johnson V (Lucey at el. 1991, hereafter L91) 
and with a recent determination of the FP in Johnson B, 
a = 1.25, (3 = 0.32, 7 = -8.895, (B98). The mean offset A7 
(defined as the mean of log Re + log /150 — a log a — /3SBe — 7 
over the data points) is 0.10 with respect to the L91 sample 
(using the slopes given in L91, i.e. a — 1.23 and /3 = 0.328) 
and 0.07 in Johnso n B. We also compar ed our data to the 
7 Samurai sample (Davies et al. 1987) obtaining A7=0.08 
and 0.07, respectively in Johnson B and V and ii) to the FP 
in Gunn r (JFK96), finding A7 ~ 0.03 with respect to the 
Coma zero point, even if the colour corrections increase in 
this case the uncertainty. The rms scatter of our points is 
0.16 in log J?e giving a standard deviation of the mean value 
of about 0.07 in A7. Considering other smaller sources of 
error such as the peculiar velocity of Coma and the colour 
transformations we can estimate the uncertainty on 7 to be 
~ 0.1. The offset found can be interpreted as a result of 
the evolution of the mean A4/L ratio, AlogA'//L ~ 0.1. In 
other words, our field galaxies are more luminous than a 
local elliptical with the same effective radius and a (dimen- 
sion and mass), consistently with the expected evolution of 



stellar populations, and in agreement with the results found 
by Bender et al. (1996), Ziegler & Bender (1997), Kelson et 
al. (1997), Schade, Barrientos & Lopez-Cruz (1997) and van 
Dokkum et al. (1998), for the cluster environment. 

If we exclude from the analysis the two galaxies E and F 
with lower signal-to-noise spectra, the mean offset becomes 
0.15 with respect to the L91 sample, 0.11 to B98, 0.13 and 
0.12 to the 7 Samurai B and V and 0.12 to JFK96. The scat- 
ter is practically unchanged (0.16). The change in the off'set 
is not negligible, but is not statistically significant. More 
data are needed to extend this preliminary study to larger 
sample in order to overcome small number fiuctuations and 
to measure the offset and scatter with smaller errors. More- 
over when a larger sample will be available, it will be possible 
to measure the offset in smaller redshift bins, thus providing 
homogeneous subsamples. 

Clusters provide a homogeneous environment for the 
formation of ellipticals. In contrast, there is no a priori rea- 
son to believe that randomly selected field galaxies have 
formed at the same epoch, and any differences in the forma- 
tion redshift will be amplified with look-back time. There- 
fore, one might have expected to find a larger scatter for our 
field relation (if any relation was to be found at all) than in 
the cluster FP. The scatter we find for our data points (0.16 
in logi?e, using the local slopes) is in fact larger than the 
scatter found by, e.g., L91 (0.075). Part of this scatter can 
be explained because we are considering a small sample of 
galaxies at different redshifts (and therefore different 7; from 
a geometric point of view, the plane is thicker because is the 
superimposition of parallel planes). It is not clear however 
whether there is a residual scatter, and a larger sample is re- 
quired to investigate the intrinsic variations of the formation 
and evolutionary history for galaxies not in rich clusters. 

6.2 Comparison with other results: The evolution 
of the Fundamental Plane 

Obviously, a larger sample of field galaxies is necessary 
to assess whether they truly obey an FP-like relation and 
whether, and how much, their properties vary with redshift. 
However, a significant body of data is now available on el- 
lipticals at intermediate redshifts, mostly from clusters, and 
it is instructive to look at our sample in the context of the 
cluster data as well. 

We consider here a total of 25 ellipticals from the pa- 
pers of van Dokkum and Franx (1996, 4 objects). Kelson 
et al. (1997, 15 objects), and from our sample (6 objects). 
The first consideration is that these galaxies, much like ours, 
fit reasonably well with the zero-redshift relations (see Fig- 
ure The zero point determined from all these galaxies 
using the L91 slopes is —8.60 (Figure ^), very similar to the 
value obtained with our six field galaxies alone. 

With this sample, it is also possible to determine the FP 
slopes a and /3 independently of the low-redshift measure- 
ments. The results depend on the method used for the min- 
imisation {a especially, while /3 is better determined). Using 
the technique described in JFK96, which equally weights 
the three FP parameters, we obtain a — 0.91 and /3 = 0.307 
(Figure We estimate the errors to be (To, « 0.15 and 
a/B ~ 0.02. This choice of the slopes reduces the scatter in 
log Re from 0.13 (with respect to the local relation) to about 
0.10. The relation a~ 10(3 + 2 = (see Section 0) still holds 
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0.91 log a + 0.307 SB,, - 7.48 




Figure 6. The intermediate-redshift (filled) points plotted 
against the L91 (Lucey et al. 1991; empty pentagons) Coma clus- 
ter sample. The filled points are from this paper (large squares), 
from van Dokkum and Franx 1996 (triangles), and from Kelson et 
al. 1997 (small squares and pentagons). All surface brightnesses 
are in Johnson V. The line corresponds to the FP relation of L91. 
The rms scatter in logiJe is 0.14. The mean offset of the field 
ellipticals is 0.10 logi?e, of the cluster sample is 0.17 




Figure 8. The intermediate-redshift points plotted against their 
own best-fitting FP relation. The points are from this paper (large 
squares), from van Dokkum and Franx 1996 (open triangles), and 
from Kelson et al. 1997 (open squares and pentagons). The rms 
scatter in log Re is 0.10. 

within the errors. If the reasoning is carried out with the 
assumption of the single parameter rj {L (x Af), this would 
lead to A»7 « —0.2. If confirmed with better statistics, this 
variation might imply a differential passive evolution as a 
function of mass, i.e. a mass-dependent formation redshift. 

In Figure ^ we show the Kormendy relation between 
SBe and Re for nine galaxies in our field sample (up to 
z— 0.647) plus the 25 other intermediate redshift cluster el- 
lipticals. The intermediate redshift data are compared in the 
figure with the L91 sample of Coma ellipticals. An average 
offset of —0.52 magnitudes is found between the interme- 
diate redshift cluster and field sample and the local Coma 
sample. The intermediate redshift galaxies are therefore on 
average brighter (at the same Re), consistently with what 
is found using the FP. In order to determine the slope ai, a 
sample with well-controlled selection biases is required. The 
intermediate redshift sample was not chosen to this aim and 
therefore we do not attempt to derive the value of ai. 

More data and carefully selected subsamples are needed 
to assess if the field Kormendy relation at intermediate red- 
shift differs significantly from the cluster one. 



Figure 7. The field intermediate-redshift galaxies parameters in 
Johnson B are shown as filled squares. The solid line is the FP 
relation given by Bender et al. 1998. The mean offset is 0.07 in 
log Re. The rms scatter in log Re is 0.15. 



7 SUMMARY 

With the instrumental setup and the procedure described 
Sections 3, 4 and 5, we have measured the Fundamen- 
tal Plane parameters of intermediate redshift field elliptical 
galaxies. As an important part of the data-reduction proce- 
dure we have carried out accurate K-corrections and spec- 
troscopic aperture corrections. Given the importance of the 
template spectrum choice for the kinematic fit, we have de- 
cided to make use of a wide set of spectra, checking how the 
results vary with the spectral type, as described in Section 
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Figure 9. Surface brightness (SBgy ) vs. effective radius (Re) for 
the same sample as in Figure 6 (filled symbols) ; the large squares 
are our data, the small squares and the pentagons are the galaxies 
from CL1358+62 and MS2053-04 (Kelson et al. 1997), and the 
triangles are the four ellipticals of CL0024 (van Dokkum & Franx 
1996). The data of a sample of Coma ellipticals (Lucey et al. 1991) 
are shown as empty pentagons for comparison, and their best- 
fit Kormendy relation is overplotted. The mean average offset is 
^ —0.5 magnitudes. Three galaxies are also shown (B, O, P) for 
which we have measured photometry and redshift but no central 
velocity dispersion. 



no significant differences in the history of the two environ- 
ments. More data are needed to perform separate fit to the 
two subsamples. 

(iii) The Kormendy relation between SBe and Re for all 
the intermediate redshift data shows an offset of ~ —0.5 
magnitudes in surface brightness with respect to the Coma 
ellipticals (L91), in agreement with what is found using the 
FP. Even though caution is needed, due to possible sam- 
ple selection biases, the variation of the Kormendy rela- 
tion suggests that the evolution found via the Fundamental 
Plane, extends to galaxies at higher redshift. Carefully se- 
lected samples of field and cluster intermediate redshift are 
needed to study the the Kormendy relation in greater detail 
and to investigate the possible differences between the two 
environments. 

Further discussion on the evolution of elliptical galaxies will 
be given in a follow-up paper (Paper II; Stiavelli et al. 1999) 
with the help of information on the line strengths. 
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4, obtaining at the same time an estimate of the systematic 
error produced by the choice of the template spectrum. 

While the number of galaxies studied here is too small 
to draw firm conclusions on the properties of field ellipticals 
at intermediate redshift, the preliminary evidence suggests 
that: 

(i) Our six field ellipticals at redshift z ~ 0.3, are in agree- 
ment with the local FP relation, with a variation of the zero 
point (A7 « 0.1), and a scatter of 0.16 in log Re - This means 
that the stellar populations of our sample of field galaxies 
are brighter than the local ellipticals, with the same size and 
mass. This data fit into a scenario in which our galaxies, at 
a look-back time of ~ 4-5 Gyrs, are evolving passively. The 
small sample and the sample selection criteria do not allow 
any more general conclusion. 

(ii) The FP obtained from o ur data and the cluster ellip- 



ticals at intermedi ate redshift (van Dokkum & Franx 1996 



Kelson et al. 1997) is well defined with a scatter of 0.13 in 
log Re with respect to the local relation. The full intermedi- 
ate redshift sample is large enough to perform the fit the FP 
coefficients independently. Using the fitting technique used 
by JFK96, we find a = 0.91 and /3 = 0.307, and the scatter 
is reduced to 0.10. By interpreting the Fundamental Plane 
as a result of homology, t he virial theorem and th e existence 
of a relation L oc M" ( Ivan Albada et al. 199 j ), different 
slopes at different redshift imply a variation of 77 with time, 
i.e., a mass-dependent evolution. The good agreement be- 
tween the field and cluster galaxies suggests that there are 
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